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Using first principles density functional theory, we study the electronic structure of multilayer 
graphene nanoribbons as a function of the ribbon width and the external electric field, applied 
perpendicular to the ribbon layers. We consider two types of edges (armchair and zigzag), each 
with two edge alignments (referred to as a- and /3-alignments). We show that, as in monolayer 
and bilayer armchair nanoribbons, multilayer armchair nanoribbons exhibit three classes of energy 
gaps which decrease with increasing width. Non-magnetic multilayer zigzag nanoribbons have band 
structures that are sensitive to the edge alignments and the number of layers, indicating different 
magnetic properties and resulting energy gaps. We find that energy gaps can be induced in ABC- 
stacked ribbons with an perpendicular external electric field, while in other stacking sequences, the 
gaps decrease or remain closed as the external electric field increases. 

PACS numbers: 71. 15. Mb, 81.05.Uw, 75.75.+a 



I. INTRODUCTION 

Recent progress in the isolation of single and multi- 
layer graphene sheets has opened up a new topic in two- 
dimensional electron systemsi"— . Electronic structures of 
single and multilayer graphene^ continue to contribute 
interesting physics and hint at possible practical applica- 
tions in carbon based electronics^"®. Recent predictions 
in multilayer graphene address chiral decomposition of 
electronic states in multilayer graphene stacks^, energy- 
gap opening only in ABC-periodic graphene stacks with 
an perpendicular external electric field applied per- 
pendicularly to the ribbon layers, and importance of 
electron-electron interactions in multilayer systems due 
to the appearance of flat bands near the Fermi level^. It 
is interesting to study whether some of these predictions 
can be extended to finite size multilayer graphene stacks 
(such as ribbons and flakes). Moreover, due to recent 
advances in fabrication of less than 10 nm wide nanorib- 
bons with the control of their edge morphologj»i2"— and 
observing magnetic edge states in few-layer graphene 
ribbon&i^ii^, our studies of interplay of magnetism and 
electric field effects in multilayer nanoribbons have im- 
portant implications in interpreting experiments. 

In this paper, we report on the electronic structure 
of multilayer graphene nanoribbons using first principles 
electronic structure methodic and address the interplay 
of magnetism, perpendicular external electric field and 
the energy gap. Our study suggests the existence of three 
classes of energy gaps in multilayer armchair nanorib- 
bons, and strong dependence of magnetic properties on 
the edge alignments and the number of layers in multi- 
layer zigzag nanoribbons. We discuss the effect of a per- 
pendicular external electric field on multilayer nanorib- 
bons and find that gaps can be enhanced in metallic 
ABC-stacked nanoribbons whereas in other nanoribbons. 



gaps decrease with increasing electric field or remain al- 
most zero. 

We begin by describing the computational method and 
the parameters used for this study in section II. In sec- 
tion III, we present the electronic structure of multi- 
layer armchair nanoribbons and discuss the width de- 
pendence of energy gaps. Band structure of multilayer 
zigzag nanoribbons and the interplay of band structure, 
edge magnetism and resulting gaps will be discussed in 
section IV. In section V, we discuss the external electric 
field effects on the energy gaps of AB- and ABC-periodic 
ribbons. In the following discussions, the direction of the 
applied external electric field is taken to be perpendicu- 
lar to the ribbon layers. Finally we present our summary 
and conclusions. 
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FIG. 1: (Color online) Schematic illustration of (a) three 
types of stacking arrangements, labeled by A, B and C, 
and (b) two types of edge alignments, Q-alignment and /3- 
alignment in multilayer graphene nanoribbons. The two edge 
alignments are distinguished by different ways of shifting the 
top layer with respect to the other in finite size multilayer 
graphene stacks. The arrows indicate the direction of edges 
along which nanoribbons span infinitely. 
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II. DENSITY FUNCTIONAL THEORY 
CALCULATIONS 

We use a plane-wave basis-based electronic struc- 
ture methodic with ultrasoft pseudopotentialsiS for core- 
valence interaction to obtain multilayer band structures. 
In our previous studies of bilayer nanoribbons^^ and 
flakesi^, it was found that the interlayer distance and 
the appearance of edge magnetism is sensitive to the 
particular local or semi-local approximation used. For 
the sake of consistency and meaningful comparisons, the 
same generalized gradient approximation (GGA)^^ was 
used to capture the edge magnetism with the fixed inter- 
layer distance of 0.335 nm. 

We note that van der Waal's interaction, which an- 
chors the layers of graphene together, is not addressed 
within Kohn-Sham DFT and recently there have been 
several attempts to include these interactions seamlessly 
in the density functionala^Si^i, which could capture the 
interlayer distance between graphene layers accurate to 
within 5%. It was shown^l, however, that although van 
der Waals interaction contribute to interlayer distance 
dependent total energy, it has weak influence on overall 
band structure at a given distance, thus we neglect the 
van der Waals interaction in this study. 

The unsaturated carbon a-orbitals were passivated 
with hydrogen atoms and the C-H distance was fixed 
at 0.1084 nm. Therefore, our results are appropriate in 
situations where ribbons are cut in a hydrogen environ- 
ment. 

In graphene sheets, there are three distinct stacking ar- 
rangements, labeled by A, B and C, as shown in Fig.[lja). 
For nanoribbons, we consider two types of edges (arm- 
chair and zigzag), each with two edge alignments namely 
a- and /3-alignments, as illustrated in Fig.[ljb). For more 
than two graphene layers, the stacking combinations mul- 
tiply and each stacking sequence presents a particular 
band structure*. Density functional theory predicts that 
periodic AB stacking (Bernal) and periodic ABC stack- 
ing (rhombohedral) are almost energetically degenerate 
in bulk graphite^. Recently, ABC stacking was realized 
on a SiC(OOOl) substrate^^ experimentally, and electric 
field-induced gap opening was predicted theoretically in 
ABC-stacked trilayer graphene as well as ABCA-stacked 
tetralayer graphene2i2ii2^. Thus, for multilayer ribbons 
with more than 3 layers, we only focus on periodic AB 
and ABC stacking sequences. 

To establish the ground state magnetic order for the 
zigzag ribbons, we tested both narrow and wide ribbons 
with non-magnetic, ferromagnetic and antiferromagnetic 
order between the layers, while we set ferromagnetic cou- 
pling along each edge and antiferromagnetic coupling be- 
tween the two edges within the same layer, as predicted 
theoretically^^. We find that interlayer antiferromag- 
netic order has lower energy than ferromagnetic, non- 
magnetic, or non-coUinear magnetic order in the ribbons. 
We also find that the same ground state is reached with 
other forms of semi-local exchange-correlation potentials 



such as PBE27, PBESoP^, RPBB^. Therefore, for cal- 
culating band structures and other related quantities in 
this article, we consider interlayer antiferromagnetic or- 
der as a magnetic ground state of multilayer graphene 
ribbons. 

The nanoribbons were placed in a supercell with vac- 
uum regions next to the width and the stacking direc- 
tion. For trilayer ribbons, we used a 1.5 nm vacuum 
region and for tetralayer, pentalayer and hexalayer rib- 
bons, 2 nm, 2.5 nm and 3 nm vacuum regions were used, 
respectively. The necessity for taking different vacuum 
sizes for increasing number of layers was guided by the 
need to avoid intercell interaction in DFT in order to 
treat the multilayer stack as an isolated system and also 
the saw-tooth type implementation of electric potential 
in the DFT code we used. 

The bulk trilayer graphene calculations^ suggested an 
external electric field, applied perpendicular to the rib- 
bon layers, as high as 4 V/nm is necessary in order to 
achieve saturation of the gap. Therefore, we use a maxi- 
mum external electric field strength of 5 V/nm to study 
its effect on the multilayer ribbon gaps. We used 68 k- 
points in the irreducible part of the Brillouin Zone (BZ) 
and kinetic energy cut-off of 475 eV. The convergence 
of the calculations was tested with respect to a denser 
k-point mesh, larger energy cut-offs as well as larger vac- 
uum sizes. 



III. ARMCHAIR GRAPHENE SHEETS 

In this section, we discuss the electronic structure of 
multilayer armchair ribbons and the width dependence 
of energy gaps. 



A. Trilayer nanoribbons 

Figure [2] shows energy band structures of armchair rib- 
bons stacked in ABA and ABC fashion. The Fermi en- 
ergy is placed at zero. Like in monolayer and bilayer 
armchair nanoribbons, both metallic and semiconducting 
trilayer ribbons are possible. Note that for the metal- 
lic ABA nanoribbon (Fig. [2][a)) we find one linear and 
one quadratic band, while for metallic ABC nanoribbon 
(Fig. [2Kc)) there is one cubic band near the Fermi en- 
ergy. The appearance of the low energy states and corre- 
sponding energy dispersions are consistent with the en- 
ergy spectrum analysis in bulk multilayer graphene^. 

We now discuss the width dependence of the ribbon 
gaps. We consider the ribbons with maximum widths 
up to 5 nm. Figure [3] shows the width dependent gap 
variation of trilayer armchair ribbons for both a- and 
/3-alignments in ABA and ABC stackings. Similarly as 
monolayer-^ and bilayer ribbonsi^, three classes of rib- 
bons are clearly seen in all cases, two semiconducting 
{N = 3p, 3p + 1) and one metallic {N = 3p -I- 2), where 
N denotes the number of carbon chains along the width 
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direction and p is an integer number. (Here, because of 
relatively smaller gap size, we call N — 3p + 2 ribbons 
as metallic though for small p, the ribbons are actually 
semi-conducting.) 

Note that compared to bilayer armchair nanoribbons, 
these gaps are consistently smaller—. As the number 
of layers increases, more energy bands appear near the 
Fermi level due to the coupling between the layers at the 
DFT level, which reduces the gap size. 



B. Multilayer nanoribbons 

For tetralayer graphene, we consider periodic AB and 
ABC stacking sequences with widths up to 5 nm. As 
with mono-, bi- and tri-layer graphene ribbons, there ex- 
ists three classes of ribbons in both stacking sequences, as 
shown in Fig. Only ribbons with a-alignment are con- 
sidered here; the /3-aligned ribbons show similar behavior 
(figures not shown). 

It is clearly seen from Fig. |4] that the gaps are consis- 
tently becoming smaller for a given width, compared to 
bilayer— and trilayer ribbons (Fig. 

The origin of this decrease in gaps may again be at- 
tributed to appearance of more energy states, near the 
Fermi level, as discussed in the previous section. Al- 
though explicitly not shown, we believe that in thicker 
ribbons (with more than four layers), the width depen- 
dence of gap may show three classes of ribbons, albeit 
with gaps smaller than the gaps in ribbons with less than 
four layers. 



an energy gap by breaking the flat band degeneracy (Fig. 

In non-magnetic /3-aligned ribbons (Fig. [SJc)), several 
flat bands appear together at the Fermi level and as a 
result, the system is again unstable and becomes mag- 
netic due to the presence of comparatively large density 
of states. When we take into account the edge magnetic 
order, an energy gap is opened at the Fermi level. (Fig. 
[5{d)). The flat bands and magnetically induced energy 
gap are also seen in bilayer graphene nanoribbons (see 
Fig. 5 and 6 in Ref. 18). The number of flat bands at 
the Fermi level or away from it depends on the number 
of graphene layers. 

We now discuss the variation of the magnetically in- 
duced gaps with the ribbon width. We consider a maxi- 
mum width of 5 nm. Figures [BJa) and (b) show a mono- 
tonic decrease of the gap with the ribbon width in both 
ABA- and ABC-stacked ribbons, respectively, each with 
a- and /3-alignments. Due to magnetism, for compara- 
ble widths, the energy gaps of zigzag ribbons are larger 
than those of corresponding armchair ribbons. For ex- 
ample, for an ABA-stacked ribbon with a width of 1.2 
nm (which translates to iV=5 and N=ll carbon chains 
for the zigzag ribbon and the armchair ribbon, respec- 
tively) and a-alignment, the gap is 0.433 eV (0.118 eV) 
for the zigzag (armchair) ribbon. The increase of the gap 
with increasing width after 3 nm in ABC-stacked ribbons 
is somewhat surprising. We increased the vacuum region 
from 1.5 nm to 3 nm and 5 nm, and repeated the calcula- 
tions for these widths but again find the same increasing 
trend. 



IV. ZIGZAG GRAPHENE SHEETS 

In this section, we study zigzag edged multilayer rib- 
bons and the induction of an energy gap due to edge 
magnetism. We first describe the interesting electronic 
structure and edge magnetism in trilayer zigzag nanorib- 
bons, and then extend our discussions to multilayer rib- 
bons. 



A. Trilayer nanoribbons 

Zigzag ribbons, due to their edge states, show magnetic 
order and their ground states are predicted to have the 
interlayer antiferromagnetic order, as discussed in Sec. 
II. Figure [5] shows the non- magnetic and magnetic band 
structure of a 1.7 nm wide trilayer zigzag ribbon stacked 
in ABA-fashion with both a- and /3-alignments. We find 
similar band structures for ABC-stacked ribbons. In non- 
magnetic ribbons with a- alignment (Fig. [Sfa)), a flat 
band appears at the Fermi level and as a result, the sys- 
tem is unstable and becomes magnetic due to the large 
density of states. Several other flat bands also appear 
away from the Fermi level. Edge magnetic order induces 



B. Multilayer nanoribbons 

We now extend our discussions to non-magnetic and 
magnetic multilayer ribbons, which show distinct fea- 
tures in their band structures depending on the edge 
alignment and whether the number of layers is odd or 
even. For illustrative purposes, we consider only AB- 
periodic ribbons with a-alignment. a-aligned ABC- 
periodic ribbons show similar behavior. For /3-aligned 
ribbons, both periodic stacking sequences show several 
flat bands appearing together at the Fermi level as in 
trilayer counterparts, and do not show any qualitative 
change with the number of layers. 

Figure [7] shows non-magnetic (left panels) and mag- 
netic (right panels) band structures of tetralayer, penta- 
layer and hexalayer a-aligned ribbons in ABA-stacking 
sequence. In non-magnetic ribbons with an even num- 
ber of layers (such as in tetralayer and hexalayer rib- 
bons), several flat bands appear away from the Fermi 
level, whereas for an odd number of layers (such as in 
pentalayer ribbons), one degenerate flat band appears at 
the Fermi level with several other flat bands away from 
it. The number of such flat bands depends on the num- 
ber of layers considered. The appearance of flat bands 
at the Fermi level results in magnetic instability and as 
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a result magnetic order induces gaps in all cases (Fig. [71 
right panels). 

With the layer antiferromagnetic order, we expected to 
find magnetically induced gaps larger in a-aligned odd- 
layered zigzag ribbons compared to the a-aligned even- 
layered ribbons due to a flat band appearing directly at 
the Fermi level. We also expected that /3-aligned rib- 
bons exhibit larger gaps compared to a-aligned ribbons 
due to several flat bands at the Fermi level. But we 
did not find any such trends. For comparison, we tested 
ribbons with layer ferromagnetic order (both a- and /3- 
alignments), and found that in /3-aligned ribbons, the 
magnetically induced gaps are consistently larger than 
their a-counterparts, albeit still no trend in gaps of odd- 
even layers is seen either in this case. 

V. ELECTRIC FIELD EFFECTS ON THE GAPS 

This section deals with the effect of external electric 
fields, applied perpendicular to the ribbons, on the con- 
finement and edge magnetism induced gaps'^°. We choose 
representative metallic and semiconducting ribbons to il- 
lustrate the effect. Figure [5] shows the variation of energy 
gap with external electric fields in ABA-stacked armchair 
ribbons with a- and /3-alignmcnts (a-b) respectively, and 
corresponding alignments in ABC-stacked ribbons (c-d). 
We consider a maximum field strength of 5 V/nm. A gap 
is induced by electric fields in ABC-stacked metallic rib- 
bons, increasing with small fields and then showing signs 
of saturation. This is consistent with the prediction of a 
gap opening in ABC-stacked bulk trilayer graphene and 
the electric field effects on the induced gaps^. In ABA- 
stacked metallic ribbons, gaps remain closed or constant 
for all electric field strengths. In all cases, for semicon- 
ducting ribbons initially with a large gap, we find mono- 
tonic decrease of gaps with increasing electric field. This 
behavior is similar to electric field effects in bilayer arm- 
chair ribbons^. 

Our calculations on bilayer nanoribbons suggested a 
critical band gap of 0.2 eV above and below which the 
external electric field change the sign of the gap values. 
Although such critical gap value is not explicitly proved 
here, we believe that a critical gap may also exist in multi- 
layer nanoribbon stacks. Therefore, we expect that these 
gaps (above the critical gap) to decrease with increasing 
strength of external electric fields applied perpendicular 
to the layers. 

To study the electric field effects in zigzag ribbons, 
we choose both wide and narrow gap ribbons. Only a- 
aligned ribbons in both stacking sequences are consid- 
ered. /3-aligncd ribbons show similar behavior (figures 
not shown). FigurelHlshows the variation of energy gap of 
the zigzag ribbon with respect to electric field strengths 



(chosen up to a maximum of 5 V/nm). In wide width 
ribbons which have initially a small gap, we find gaps in- 
creasing with increasing electric field strengths for ABC 
stacking while for ABA stacking, gaps remain closed or 
constant. In narrow width ribbons which have initially 
a large gap, the gap decreases with field. This trend is 
similar to that observed in bilayer zigzag nanoribbonai^. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have studied the electronic proper- 
ties of armchair and zigzag multilayer graphene nanorib- 
bons both with and without external electric fields using 
first principles density functional-based electronic struc- 
ture method. We consider both AB-periodic and ABC- 
periodic nanoribbons with two different edge alignments, 
referred to as a- and /3-alignments. Armchair edged mul- 
tilayer ribbons exhibit three classes, two semiconduct- 
ing and one metallic. The energy dispersion of metal- 
lic armchair nanoribbons near the Fermi energy is con- 
sistent with the chiral decomposition in bulk multilayer 
graphene. The gap in multilayer armchair ribbons, for a 
particular width, is found to be smaller than in the cor- 
responding bilayer ribbons. We find that magnetically 
induced gaps in multilayer zigzag nanoribbons and the 
modulation of the gap values depending upon the type of 
edge alignments. ABC-periodic nanoribbons with a small 
gap show enhancement in the gap values with increas- 
ing electric field strengths, while AB-periodic nanorib- 
bons with a small gap show no increase in the gap values 
with increasing electric fields. For nanoribbons with a 
large gap, the gap values decrease with increasing electric 
fields. In view of encouraging advances in the fabrication, 
control of graphene edges and observation of magnetic 
edge states, we believe that our studies of magnetism 
and electric field effects will be important for designing 
graphene based nanodevices. 
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FIG. 2: Energy band structure of trilayer armchair nanorib- 
bons for (a) ABA-a, metallic nanoribbon, (b) ABA-a, semi- 
conducting nanoribbon (c) ABC-a, metallic nanoribbon, and 
(d) ABC-Q, semiconducting nanoribbon. The chosen widths 
for metallic and semiconducting ribbons are 0.86 nm and 1.11 
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FIG. 3: (Color online) Variation of the energy gap with 
widths of (a) ABA-a, (b) ABA-/3, (c) ABC-q, and (d) ABC- 
/3-aIigned nanoribbons. Three classes of the nanoribbons are 
clearly seen in (a)-(d) and specified by A'' = 3p, 3p + 1 and 
3p+2 where A is the number of carbon chains along the width 
direction. Here p =1,2,- • • ,13, which translate to nanoribbons 
with widths up to 5 nm. 
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FIG. 4: (Color online) Variation of the energy gap with 
widths of tetralayer nanoribbons (a) ABAB-q, (b) ABCA-a. 
Three classes of nanoribbons are clearly seen. These classes 
are specified by = 3p, 3p + 1 and 3p + 2 where N is the 
number of carbon chains along the width direction. Here p 
— 1,2,- ■ ■ ,13, which translate to nanoribbons with widths up 
to 5 nm. 
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FIG. 6: Variation of energy gap with the widths of (a) ABA 
and (b) ABC-stacked trilayer zigzag nanoribbons. 
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FIG. 7: Energy band structure of multilayer zigzag nanorib- 
bons for ABA-stacked (a) tetralayer (b) pentalayer and 
(c) hexalayer with a-alignmcnt. The left panels show the 
nanoribbons with non-magnetic edge atoms and the right pan- 
els show the nanoribbons with edge atoms ordered antiferro- 
magnetically layerwise. The width is chosen as 1.7 nm. 



12 



0.8 
0.4 


0.8 

0.4 



o. 
O 



0.6 



0.3 




0.9 



0.6 



0.3 



(a) ABA-a 



O W=l . 1 1 nm 
aW=0.8 6 nm 



(b) ABA-p 



O W=l . 1 1 nm 
aW=0.86 nm 



(c) ABC -a 



O W=l . 1 1 nm 
aW=0.8 6 nm- 



"E-- .,^ A -- - -a 

...A 13 _ ^ i 



"■--o-----e — e--'-'-~-if) 



(d) ABC-P 



O W=l . 1 1 nm 
aW=0.8 6 nm 



1 2 3 4 5 
Eext(V/nm) 

FIG. 8: (Color online) Variation of the energy gap with per- 
pendicular external electric field for ABA-stacked and ABC- 
stacked trilayer armchair ribbons. Panels (a) and (b) are for 
a- and /3-aligned ABA-stacked ribbon respectively whereas 
panels (c) and (d) are, respectively for ABC-stacked a- and 
/3-aligned ribbons. Both metallic (W = 0.86 nm) and semi- 
conducting (W=l.ll nm) nanoribbons are considered which 
corresponds to A'^=8 and 7V=10 respectively, where A'^ is the 
number of armchair carbon chains along the width direction. 
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FIG. 9: (Color online) Variation of the energy gap with per- 
pendicular external electric field for (a) ABA and (b) ABC- 
stacked trilayer zigzag nanoribbons with a-alignments. /?- 
aligned nanoribbons show similar behavior. Both narrow (W 
= 1.74 nm and 1.85 nm) and wide (W= 5.11 nm and 3.75 
nm) nanoribbons are considered. 



